www.nature.com/scientificreports 


scientific reports 


OPEN 


| ® Check for updates 


The computational model 

of nanofluid considering heat 
transfer and entropy generation 
across a curved and flat surface 


Sayer Obaid Alharbi?, Florentin Smarandache?, Awatif M. A. Elsiddieg?, 
Aisha M. Algahtani*~, M. Riaz Khan®, V. Puneeth® & Nidhal Becheikh” 


The entropy generation analysis for the nanofluid flowing over a stretching/shrinking curved region is 
performed in the existence of the cross-diffusion effect. The surface is also subjected to second-order 
velocity slip under the effect of mixed convection. The Joule heating that contributes significantly to 
the heat transfer properties of nanofluid is incorporated along with the heat source/sink. Furthermore, 
the flow is assumed to be governed by an exterior magnetic field that aids in gaining control over 

the flow speed. With these frameworks, the mathematical model that describes the flow with such 
characteristics and assumptions is framed using partial differential equations (PDEs). The bvp4c solver 
is used to numerically solve the system of non-linear ordinary differential equations (ODEs) that are 
created from these equations. The solutions of obtained through this technique are verified with the 
available articles and the comparison is tabulated. Meanwhile, the interpretation of the results of 

this study is delivered through graphs. The findings showed that the Bejan number was decreased 

by increasing Brinkman number values whereas it enhanced the entropy generation. Also, as the 
curvature parameter goes higher, the speed of the nanofluid flow diminishes. Furthermore, the 
increase in the Soret and Dufour effects have enhanced the thermal conduction and the mass transfer 
of the nanofluid. 


List of symbols 

Re; Reynolds number 

bh Dynamic viscosity [kg m~!s~}] 

7) Nanoparticle’s concentration 

v Kinematic viscosity [m/s] 

s Arc length coordinate with respect to the curved surface 
(pCp) Heat capacity 

A Buoyancy ratio parameter 

a Thermal diffusivity [m?s~!] 

Cf Coefficient of skin friction 

Shy Sherwood Number 

Nu; Nusselt Number 

Pp Dimensional pressure [kg m~!s~7] 

fa Temperature of the fluid [K] 

K Thermal conductivity [Wm IkK-} 

a Constant related to stretching and shrinking of the sheet 
p Density [kgm >] 
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Trs Wall shear stress 

Of Os Electrical conductivity [Sm7]] 

A Mixed convection parameter 

S Suction parameter 

K Dimensionless curvature parameter 

u,v Velocity components in s, r directions respectively [ms~]] 
R Curvature of the curves belt 

Cc Concentration of the fluid [molm~7] 

Cy,Coo Concentration near and far away from the surface respectively 
Ty, Too Surface temperature and temperature far away from the surface respectively 
r Normal to the tangent at any point of the curved surface 
M Magnetic field parameter 

Bo Magnetic field strength [T] 

Dm Molecular diffusivity [m2s~]] 

Pr Prandtl number 

jw Wall heat flux 

ky Boltzmann constant [8.314 J/Kmol] 

Vy Suction velocity 

n Fitted rate constant [W m~?K~4] 

k* Mean absorption coefficient [m7] 

P Dimensionless pressure 

Ec Eckert number 

Cs Concentration susceptibility [kgm~ >] 

dw Wall heat flux 

ky Chemical reaction rate [mol L~!s~?] 


* 


Oo Stefan-Boltzmann constant [Wm7?K~*] 
Eq Activation energy [J] 

kr Thermal diffusion ratio 

f Stream function 

n Similarity variable 

f' Dimensionless velocity 

6 Dimensionless temperature of fluid 

0) Dimensionless concentration of fluid 

L Slip length 

Br Brinkman number 

H Diffusion parameter 

Rd Radiation parameter 

Du Dufour number 

Sc Schmidt number 

E, Dimensionless activation energy parameter 
@ Temperature difference parameter 

€ Dimensionless slip length 

T Dimensionless chemical reaction rate parameter 
QO, Concentration difference 

E1; €2 Velocity slip parameter 

Sr Soret number 

Ne Local entropy generation 

y Stretching/shrinking parameter 


The term nanofluid was coined in the year 1995 by Choi’ considering the fact of suspending the nanoparticles 
would enhance the capacity of the conduction of heat of the regular fluid. Nanofluid is a type of heat carrier 
that consists of a base fluid and metal particles of size 10~? m. The nanoparticles minimal surface area allows 
the enrichment of heat/energy transfer. It is widely employed in various fields like automobiles, nuclear reac- 
tors, refrigerators, cooling of electronic appliances, and many other household essentials due to the nanofluid’s 
increased thermal conductivity. Also, these nanofluids find application in various manufacturing industries. In 
this regard, Khan and Puneeth’ investigated how Brownian motion and thermophoresis affected the thermal 
characteristics of nanofluid. Sharma et al.’ gave a piece of collective information on the recent advances in 
machine learning that helps in utilizing it in the analysis of thermal properties of nanofluids. Zhang et al.* ana- 
lyzed the bioconvection process in enhancing the nanoparticle distribution in the nanofluid. Pramuanjaroenkij 
et al.° completed a numerical study to understand the behavior of various thermal conductivity models for fluid. 
Further, Puneeth et al.° came to the conclusion that as the Casson parameter is increased, jet speed of the Cas- 
son nanofluid decreases. A comparative analysis was carried out by Bheshti et al.’ for the flow of a nanofluid in 
an annulus. Algahtani et al.* observed an enhancement in the temperature profile of nanofluid flowing across 
a cylinder for higher radiation under the action of viscous dissipation. The flow of sodium alginate suspended 
with Al,O; and Cu was studied by Nadeem et al.’ using the fuzzy hybrid nanofluid model. Atashafrooz'®" 
studied the dynamics of water suspended with nanoparticles to understand the effectiveness of water with solid 
suspensions as the heat carrier. Kumar et al.” performed irreversibility analysis of an unsteady non-Newtonian 
Micropolar fluid containing CNT to analyse its thermal features. Maiti et al.'° implemented the fractional order 
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model to study the heat transfer properties of blood under the influence of thermochemical effects. Dhlamini 
et al.'* deliberated the phenomena of bioconvection in the flow of nanofluid past a hot surface. Atashafrooz 
et al.'° employed simulation to realize the pattern of the nanofluid flow in the interior of a trapezoidal enclosure. 
Furthermore, they’’ considered the impression of Lorentz force on the velocity of the nanofluid flow using math- 
ematical model and conclude that the increasing strength of Lorentz force decreases the velocity. More studies 
related to nanofluid can be read in!”~”. 

The level of irreversibility occurring during a process is determined by entropy generation which is described 
using the second law of thermodynamics. This law helps in minimizing the entropy generation that enables the 
identification of optimal engineering system designs. Meanwhile, entropy generation can be used as a criterion 
to examine the working of engineering appliances. Due to these advanced features of entropy generation, many 
researchers have theoretically analyzed its effect on the heat transport properties of nanofluid. For instance, Alsu- 
lami et al.” studied the heat transfer in the swirling flow of nanofluid using modified Kriegger-Dougherty model. 
Sarada et al.*’ analysed the impact of exponential heat source on the thermal properties of ternary nanofluid. 
Zhang et al.”” used the Joule heating to examine the impact of the magnetic field on the production of entropy 
in nanofluid. Khan et al.”’ discussed the thermal features of a Casson nanofluid flowing across an expanding 
sheet. Alsulami et al.** analysed the non-equilibrium conditions for the nanofluid flow comprised of Tig Al, V and 
AA7075 nanoparticles. Punith et al.’ studied the impact fo the induced magnetic field that is generated due to the 
flow current. Huang et al.”° analysis the friction drag caused by the Lorentz effect on the flow of nanofluid across 
a curved surface. Shoaib et al.’” discussed the process of entropy generation in detail for the flow of nanofluid 
and hybrid nanofluid respectively across a stretching surface and rotating system. 

The inclusion of nanoparticles into the fluid not only enhances its heat conduction capacity but also influ- 
ences the fluid in many different aspects. For instance, the viscosity of the nanoparticles significantly increases 
based on its physical properties, density is often increased and as discussed earlier, the thermal conductance 
also increases. The difference in the thermal expansion coefficients between the nanoparticles and base fluid 
will cause thermal expansion induced convection which in turn affects the heat transfer. Jamshed et al.”* studied 
the entropy effect on the flow of second grade nanofluid. Atashafrooz et al.” studied the entropy generation 
and the impact of Bejan number along with the analysis of the thermal features of nanofluid. Mandal et al.*° 
considered the features of heat transport of a nanofluid considering the Entropy generation with the existence 
of microorganisms. Oyelakin et al.*! obtained an optimized model for accurate estimation of entropy generation 
for the Casson nanofluid. Nayak et al.*’ designed a 3D model to analyze the impact of radiation and the Entropy 
generation over the flow features and heat transfer of a nanofluid. 

The addition of a temperature gradient at the borders causes a more dynamic and significant effect to be 
produced in the fluid flow. This non-homogeneity in the thermal distribution produces a buoyancy effect which 
alternatively will have an impact on the coupled fields of velocity and temperature in the medium. The knowledge 
of mixed, forced, and natural convective flow plays a very important role in the fluid dynamics point of view as 
well as in practical engineering applications. The mixed convection impact on the flow of hybrid nanofluid was 
analyzed by Xia et al.*° with multiple slips at the boundary. Dawar and Acharya™ studied the timed dependent 
flow of nanofluid influenced by a mixed convection. Wang and Xu” obtained a very accurate solution for ana- 
lyzing the influence of convection over the nanofluid flowing in a lid-driven cavity using the wavelet-homotopy 
method. Khan et al.*° demonstrated that as the mixed convection parameter increases, the velocity of the hybrid 
nanofluid decreases. Tian et al.*” showed that mixed convection impacts positively on the nanofluid flow in an 
inclined cavity. Wahid et al.** concluded that the strength in the magnetic field diminished the velocity of the 
flow across a porous vertical cone. It was observed by Ketchate et al.*’ that the nanoparticles of blade shape 
stabilized the convective flow. Mahmood et al.*° designed a mathematical model to interpret the heat flow and 
the motion of nanofluid in a square cavity. Further Muhammad et al.*!~® elaborated the use of openFOAM in 
analyzing the fluid properties. 

The nanofluid flowing across a stretched sheet has a significant role in practical applications including the 
manufacturing of glass fiber, plastic film extraction, condensation of liquid films, paper production, etc. In these 
applications, a large amount of heating is involved and hence cooling of appliances becomes a mandatory process 
to maintain an optimum temperature. Thus, many scholars are actively working on analyzing the thermal/energy 
characteristics of the nanofluid motion across the stretching sheet and other various geometries. For instance, 
Reddy et al.** studied the significance of radiation on the stagnation point flow of nanofluid over a curved 
surface. Abbas et al.*° completed the numerical investigation to estimate the heat transfer of time-dependent/ 
unsteady flow of micropolar fluid flowing past a curved region. Qian et al.*° framed a mathematical model that 
described the significant role of the Lorentz force caused by the magnetic field on the micropolar fluid flowing 
across a curved stretching sheet. Naveen et al.*” incorporated the model of Cattaneo-Christov to design the heat 
flux in the flow of nanofluid going through a curved stretching sheet. Khan et al.** considered the impact of gold 
nanoparticles on enhancing the thermal properties of blood flowing through a curved surface. Ashraf et al.” 
examined the enhancement in the heat transfer rate of a fluid suspended with Al,O; and Fe;0, nanoparticles. 
Hayat et al.°° deliberate the entropy generation process in the nanofluid going through a curved region. To ana- 
lyze the heat and mass transport of a nanofluid passing through a curved sheet, Imtiaz et al.*! created the energy 
equation using the Soret-Dufour model. Alblawi et al.** applied the Buongiorno’s model to analyze the effect of 
the major slip mechanisms on the heat transport features of nanofluid. 

The detailed literature review provided above indicated that the availability of resources on the flow of nano- 
fluid past a curved stretching/shrinking surface is limited. Meanwhile, there were no resources available which 
can describe the influence of activation energy over the motion of nanofluid across a curved surface by consider- 
ing the Soret and Dufour effect. Thus, the authors have incorporated these effects along with the Joule heating, 
second order velocity slip, heat source/sink radiation and the magnetic field effects. The practical applications of 
the Soret effects includes the isotope separation, purification of gases and liquids, isoelectric focusing, purification 
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and analysis of proteins, biomolecules, and drug molecules, whereas the applications of Soret effect includes heat 
exchangers, combustion processes, cryogenic systems, and semiconductor processing. The mathematical model 
for the considered effects is constructed using PDEs, and this system of equations is subsequently translated to 
yield the proper system of ODEs. The solutions to the resulting system of ODEs are obtained using the MATLAB 
bvp4c tool, and the outcomes are displayed graphically using graphs. 


Basic governing equations 

Consider the dissipative mixed convective flow of a nanofluid across a stretching/shrinking curved sheet as 
shown in Fig. 1. The flow is two-dimensional and the fluid is incompressible including the effect of Joule heating, 
thermal radiation, second order velocity slip, activation energy and heat generation/absorption. The Dufour and 
Soret numbers were also appropriately considered in the energy and mass diffusion equations. The two direc- 
tions, r and s were taken as being respectively vertical to the surface and along the surface with the stretching/ 


shrinking and second order slip velocity u = as + Ly (# 7 ) + Lo (3 +, —1 ou ) as well as the 


u 
r+R ar? r+R Or (r+R)? 
free stream velocity is u — 0. Note that a = 0,a < 0 anda > 0 correspondingly indicates the static, shrinking 
and stretching surface, whereas Lj and L> respectively signifies the first and the second order slip coefficients. 


There was a fixed radial magnetic field with an intensity of Bp. Given these factors, the following are the govern- 


ing boundary layer equations’. 

or + Rv] + RH =0 (1) 
r+R)v]4 = 0, 

or ds 
10 5 
—P__* 20, (2) 
por r+R 

R Ps (vs + Ru len ) & 1 du u ) 
v = 
r+Rdas  \"ar r4+R0s r4+R V\ or r4R Or (r +R) (3) 


+ gpBe(T — Too) + geBr(C — Coo) — o Bou, 


dT = Ru =) k (#T 1 oT o2, 1 (eT 1 2) = 
v = a Bovu' 4 
dr r+R 9s pCp \ Or? or +R Or pCp pCp 3k* 


Hu (# u y+ Q P+ 1 Puke (4 1 ac) 
T Co T > 


dr2 or +R Or 
pCyp\ dr r+R pCp pCy Cs dr2 oo r+RoOr 
(4) 


ac Ru 0c aC 1 aC Dincskr ( 0?T 1 OT K2 ( T \" _ ta 
=D i! RT (C — Co). 
Gas. =) (et eR x) Lon Gatrem) (ge) : ( ve) 


(5) 


The related boundary conditions are specified as 


(1, v) 
u-0 T -T,,C -~C. 


Nanofluid 


v=—Vy du Pu 1 du u 
Or? r+ROr (r+R)? 


O 


Figure 1. Schematic diagram of the problem displaying two-dimensional motion of the nanofluid across a 
curved stretching/shrinking surface. 
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(6) 


_ a 1a = 
u=as+1i(# - 4) +o(& | FER Or 5). Yw> 
T=Ty,C = Cy,atr =0 

u > 0, 94 > 0,T > To,C > Cogasr > 00 


The following dimensionless transformations can be utilized to transform Eqs. (1-6) into a non-dimensional 


structure. 
u = bsf'(n)in = \/2rv = — Be Vbuf (n).p = pb? P(n) \ 6 
= 6(n)(Tw — Too) + Too, C = (0) (Cw — Coo) + Coo 

Asa result, the following structure is assumed by the resulting non-dimensional equation. 
P11 in 
an 4K?” 8) 

2 K 2 K 
p=f"” " / = " 464 A* M2 
7+K)S aoe al yeRY? n+Ke 7 aoe” emer: 


(9) 
! K , Td 7 7 2 y\2 (0 1 ‘ 
6 ——fo'+E,| f” — —— ME, 6 +D ——4’'} =0, 
)+ + ( 7.) + (f')" + QO + ule" + ee 
(10) 
o” + 4 +54 fo! + Sr 6” 4 9! — t(1 + 6)"Exp| — - =0, (11) 
n+K n+K n+kK 1+ 6 


We may remove the pressure P from Eq. (9) based on Eq. (8). Therefore, the sum of Eqs. (8) and (9) can be 
expressed as. 


1 
—(1+Rd){ 6” + —— 
a (l+ Ray (0" + 


iv 2 A 1 / K wt i gl aA 12, 
+ t + 
Fe dake mee ae KO if) ae! f) fies 
=4 moet es: +6) ou +o’) -ae(p+ ne) =" 
U] 
By re-arranging Eqs. (10) and (11) we get 
1 1 1 1 K 
—(1 + Rd) — SrScDu $0" —(1+Rd D 6’ + QO 
in + Rd) — SrSc ub +t a + Pe SrSc. ae i} Q 
K ri ) 
DuS ‘LE w 13 
+ te ae * (yf 7+K (13) 


2 Ey 
+ M°E,(f’)” + tDu(1 + wo)" ScExp( — 7 =) =0, 


$” + — <9! + sce] * fo! +. sr(0” + 6") — 200 4 0)"Exp( ——*! =0, (14) 
n+K ' \q+K "| n+K pues . 


The boundary conditions defined in (6) are non-dimensionalised to the following form: 
fO=SFO = tar{f"O+LO} +f" +42 - Oh, 


9(n) = 1,¢@) =1 ; (15) 
f'n) = 0 f'n) = 0,8(n) = 0,6(m) = Oasn > 00. 


Below are definitions for the dimensionless parameters that result from Eqs. (12-15). 


= a 2_ oBy? 3 _ G 79x _ G _ pvp = uw 

K=R,/2,M ap 0A = peed” = ga Pr po Ee = OT=Tay 
_ 160*T3 py, — Dink r(Cw—Coo) a Dincsk(Tw—Too) ., a 

Rd = ape, Du = CoG) Ty Tas)? Sc= 5 Sr= Tmv(Cw— - Y= ¢-: (16) 
— KF Tw—Teo = 

T= Cl po= “Ee By = 56) = Ly y/5,62 =1n4,Q= ae 


Entropy generation modeling 
The entropy generation is defined as 


é k 1 160*T 50 aT\? a fau  u \?  oBR, , RDm / IC * RD aT aC 
gen TT kk Or) ” To\ Or” Ror) “Be” * Ge Vor Tx \ar ar)’ 


(17) 
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* 73 « 2 
In above equation, A {1 + I ee \ () represents the heat transfer irreversibility, 
2 foe) 


LL ou u < ais 3 : A Pee 
ls or + Ray) Tepresents the viscous dissipation irreversibility, 
20 17 represents the Joule heating irreversibility, 

T P 8 
RDin (3c)? 4. RDm ( aT ac 
Coo \ Or Too dr or 

eep in min at Ris the universal gas constant. 
Keep d that R is th lg tant 


Applying Eq. (7) results in the dimensionless version of Eq. (17), which may be expressed as 


represents the mass transfer irreversibility. 


4 1 ? 
Ne= (1 aa 5Rd) wo!” + B, (r" at =") MBF? + H— 9" + H6'¢’, (18) 


where the definitions of Ng, B,, H, w and w are given below. 


TooVSG pbs? RDm(Cw — Coo) Tw—To AT Cyw— Cos — AC 
NG = ry r = > => > 2) —4 = > (ony = = > 
bkAT kAT k Too Too Coo Coo 
(19) 
The definition of the dimensionless Bejan number formula is 
__ Entropy generation associated to heat and mass transfer 
— Total entropy generation , 20) 
This implies that. 
ae (1+ $Rd)w0” +26? + HO'6' on 


5 
(1+ 4Rd) wo +B,(f" + okt’) + MB,f’? + H2.6” + H0'9’ 


Solution method 

The above dimensionless Eqs. (12-15) form boundary value problem and hence these equations are changed 
to a system of equations such that it forms an initial value problem. This enables us to implement the bvp4c 
numerical method using MATLAB where the interval of integration is assumed to be in the range 0-5 about 
the mesh point 70. Further, the accuracy is set to 10~® and the following relations are described to perform the 
conversion to initial value problem: 


fin) = yO) f'n) = yQ2)f") = yO)F'"(n) = WA) = YW 
A(n) = y(5),6'(n) = y(6),4"(n) = yy 


$(n) = y(7), ¢'(n) = (8), 6" (n) = yy3 


we could rewrite the resulting Eqs. (12)-(15) as 


2 1 1 
= Fi 3 2 Dy(4) — yQyB 
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_3f 72) «f 97) ie 1 
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oo RD sa | 4K ae nik’! i} 
K ] 2 
y(6) — Qy(6) — DuSe—yDy(8) — Fe (1 - ae +p) (23) 


aie 2 = n Fi 
M E.(y(2)) tDu(1 + wy(5)) scBsp( =5)} 


K 1 r E 
I y(8) se{ = pK D8) sr(y, + 7) —t(1+ @y(5)) exp(- >) \ 


n+kK 1+ wy(5) 
(24) 


The boundary conditions corresponding to (15) will take the following form: 
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Saba et al.°° | Ishaket al.°” | Mishra et al.°* | Present study 

Pr —6'(0) 

0.72 0.80884 0.8086 0.8088 0.808688 

1 1.00001 1.0000 1.0000 1.000018 

3 1.92368 1.9237 1.9236 1.923599 

7 3.07226 3.0723 3.0723 3.072316 

10 3.72068 3.7207 3.7206 3.720589 

100 12.29407 12.2941 12.2941 12.294108 


Table 1. Numerical outcomes of —6’(0) against diverse values of Pr providing the validity of the current work. 
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Figure 3. Influence of ¢2 0n the f’(7). 
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Figure 4. Influence of K on the f’(7). 
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Figure 5. Influence of / on the f’(7). 
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Figure 6. Influence of 2* on the f’(7). 
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Figure 7. Influence of Du on the 6(7). 
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Figure 12. Influence of Sc on the #(7). 
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Figure 14. Influence of Br on the Ng. 
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Figure 16. Influence of ¢2 on the Ng. 
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Results and discussion 

The mathematical model framed using the PDEs as specified in (1)-(6) were transformed to ODEs (12)-(15) 
using the transformation given in (7). The subsequent system of transformed ODEs was solved by imple- 
menting the bvp4c package as described above with an accuracy of 10~° and the solutions are validated by 
comparing with the existing literatures. The analysis is performed to understand the significant role of dif- 
ferent features of fluid flow parameters on the mass and heat transfer profiles of the nanofluid. Table 1 
and the graphical results (Figs. 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21) are used to 
report the model’s outcomes. The graphs were obtained by varying one or two of the fluid parameters 
and keeping the rest as constants. The constant values of the parameters for this study were chosen to be: 
n=1,S=2,M =0.1,R =0.1,Pr = 6.2, Ec = 0.1, Du = 0.2, Sc = 2,t = 0.1, £, = 04,4 = 0.1, 7* = 0.1, 

€) = 0.1,€2 = 0.01, K = 100,Sr = 0.4,Q=1. 

The effect of the both the slip parameters on the flow rate of the nanofluid along a stretching surface is 
depicted to be decreasing in the Figs. 2 and 3 respectively. The increase in these parameters creates a back flow 
at the boundary region which opposes the fluid. Moreover, ¢1 has opposite effect across a shrinking surface 
whereas €2 has same decreasing effect. Thus, the decrease is observed for higher values of slip parameters €2 across 
a stretching and shrinking surface whereas €, possess opposite effect across a stretching and shrinking surface. 
Figure 4 illustrates the direct relationship between the curvature parameter K and the sheet radius. As the sheet 
radius lowers, less area is available for particles to adhere to, which in turn causes the stretching and shrinking 
rate to reduce and the fluid velocity to decrease. The higher values of mixed convection parameter signify a greater 
temperature gradient indicating a lighter density of fluid. Thus, it enables the fluid to flow at a faster velocity 
across a stretching surface as shown in Fig. 5. Although opposite effects were shown across a shrinking surface. 
The higher buoyancy ratio parameter reduces the temperature gradient and the concentration gradient increases. 
This makes the fluid flow difficult and hence a diminishing effect is observed in Fig. 6 for the increase in 1*. 

As the higher Du, there is an increase in temperature and thermal diffusion, as seen by the accelerating 
behavior of 6 (7) in Fig. 7. The heat generated due to the internal friction occurring because of the flow boost with 
the rise in the Eckert number. The boost in the radiation signifies that the heat dissipated through the surface is 
more which is conducted by the nanofluid hence the amount of heat that a nanofluid conducts will be more as 
shown in Fig. 8. The growth in the suction parameter diminishes the fluid temperature as shown in Fig. 9. The 
upsurge in the Soret number signifies a greater temperature difference which shows that the temperature at the 
surface is evidently high than the surrounding. The temperature of the nanofluid upsurges overall as a result of 
the nanofluid absorbing this high heat, as depicted in Fig. 10. 

As Sr increases, the concentration profile ¢(7) also increases, as seen in Fig. 11. The "effect of ratio of tem- 
perature difference to concentration difference" is what Soret number is defined as. This makes it clear that 
a larger concentration profile #(7) is produced by diffusive species with higher Soret values. The increasing 
Schmidt number indicates that the concentration diffusion is more as compared to the previous stage, and this 
results in the decrease in the nanoparticle concentration profile as shown in Fig. 12. Enhancing the values of 
non-dimensional activation energy E; describes the cumulative response in the @(7) as shown in Fig. 13. The 
mathematical relationship in Eq. (16) clarifies that the low temperature and high activation energy decreases the 
rate of chemical reaction, which trigger the chemical reaction process to slow down. Asa result, the concentra- 
tion @(n) of the nanofluid enhances. 

The ratio of the viscous dissipation to the external heating is termed Brinkman number and the increase in 
this value enhances the viscous dissipation and slows down the conduction of heat produced. Thus, decreasing 
the entropy generation as shown in Fig. 14. Similarly, the rise in the velocity slip parameters ¢) and €2 enhances 
the local entropy generation as shown in Figs. 15 and 16 respectively. 
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Figure 17 indicated that the rise in the parameter concerning the curvature of the radius shall diminish the 
entropy profile because of the change in the boundary’s physical configuration in the radial direction. The increas- 
ing values of w signifies a larger difference in the concentration of nanoparticle which reduces the entropy profile 
as exposed in Fig. 18. With the rise in Brinkman parameter, the pressure drop also diminishes and as a conse- 
quence, a reduction in the Bejan number is recorded as shown in Fig. 19. Whereas the velocity slip parameters 
€,and €2 increase the Bejan number as shown in Figs. 20 and 21 respectively. Meanwhile, the cumulative values 
of the chemical reaction parameter reduces the nanoparticle concentration which further reduces the pressure 
drop and as a consequence, the Bejan number declines effectively as seen in Fig. 22. 


Conclusion 

The RKF-45 method is used to numerically analyze the Entropy generation for the flow of nanofluid across a 
curved stretching/shrinking surface. The mathematical model was framed by assuming the flow to be steady 
and laminar subjected to Magnetic field and mixed convection. Furhtermore the enerfy equation was designed 
using the cross diffusion model and the mass transfer equation is equipped with the Arrhenius activation energy. 
With these factors, the system of PDEs that governed the flow was transformed into a system of ODEs and the 
subsequent system of equations were solved by implementing RKF-45 method. The conclusions of the study 
indicate that: 


e With increasing velocity slip parameters, there was a drop in the flow velocity, the Bejan number, and the 
local entropy generation. 

e The higher values of the mixed convection parameter respectively increased and decreased the speed of 
the flow across a stretching and shrinking surface whereas the higher values of buoyancy ration parameter 
decreased the flow speed. 

e A direct relationship is present between the fluid velocity and the curvature parameter. 

e The Soret effects contributed in enhancing the total heat conduction by the nanofluid while the greater values 
of Eckert number showed a decreament in the thermal conduction. 

e The higher Schmidt numbers decreased the mass transfer profile whereas the higher Soret numbers enhanced 
the concentration in the nanofluid. 

e ‘The increase in the Brinkman number showed increasing effects on the local entropy generation whereas it 
had a diminishing impact on the Bejan number. 

e A reduction in velocity is detected for the rising values of slip parameter across a stretching and shrinking 
surface. 

e The mathematical model presented in this article is completely based on the Nvier stokes equation which 
often assumes ideal conditions and might not reflect the real world situations. Therefore, this model can be 
adapted to understand the flow complexity and make necessary modifications as per the requirements. 

e Enhancing the values of E describes the cumulative response in the (7). 


Data availability 
Upon reasonable request, the datasets utilized and/or analyzed in the current work will be made available by 
the corresponding author. 
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